Tensor polarizability and dispersive quantum measurement of multilevel atoms 
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Optimally extracting information from measurements performed on a physical system requires 
an accurate model of the measurement interaction. Continuously probing the collective spin of an 
Alkali atom cloud via its interaction with an off-resonant optical probe is an important example 
of such a measurement where realistic modeling at the quantum level is possible using standard 
techniques from atomic physics. Typically, however, tutorial descriptions of this technique have 
neglected the multilevel structure of realistic atoms for the sake of simplification. In this paper 
we account for the full multilevel structure of Alkali atoms and derive the irreducible form of the 
polarizability Hamiltonian describing a typical dispersive quantum measurement. For a specific set 
of parameters, we then show that semiclassical predictions of the theory are consistent with our 
experimental observations of polarization scattering by a polarized cloud of laser-cooled Cesium 
atoms. We also derive the signal-to-noise ratio under a single measurement trial and use this to 
predict the rate of spin-squeezing with multilevel Alkali atoms for arbitrary detuning of the probe 
beam. 

PACS numbers: 03.65.Ta, 42.50.Lc, 02.30.Yy 
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I. INTRODUCTION 



Information gained by performing measurements on a 
quantum system can reduce uncertainty about one or 
more of its physical observables. It is, however, a basic 
property of quantum mechanics that measurements are 
invasive in the sense that they necessarily degrade one's 
ability to make subsequent predictions about the values 
of complementary observables 0|. This type of distur- 
bance is often called measurement backaction, and it is 
a natural consequence of the Hamiltonian coupling be- 
tween a probe (such as an electromagnetic field mode) 
and the system of interest. In a special class of "back- 
action evading" experimental scenarios, it is possible to 
channel the disturbance into observables that are not 
dynamically coupled to the main quantities of interest. 
When such measurements are performed with minimal 
technical imperfection on systems whose initial prepara- 
tions are sufficiently pure, which qualifies them as what is 
referred to in the quantum optics literature as quantum 
non-demolition (QND) measurement 0,01 it is possible 
to create conditionally squeezed states of the measured 
observable. 

While measurement-induced squeezing can easily be 
understood in an abstract sense, predicting the precise 
degree of squeezing that can be achieved in a realistic 
experiment requires detailed physical modeling of the 
system-probe interaction (in addition to any operative 
decoherence mechanisms). Squeezed states of atomic 
spins have recently emerged 0, H, Ifl H a versatile and 
robust resource for quantum information science 
and quantum metrology |0, ^> Qi Ql- In these 
contexts, the degree of spin squeezing is directly linked to 



entanglement measures, to achievable reductions in aver- 
aging times for precision measurement, and to achievable 
improvements over communication protocols that utilize 
only classical information resources. 

Theoretical analyses of measurement-induced spin 
squeezing typically consider a system ol N ^ 1 atoms 
whose collective spin is described by an observable 
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where — 
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(1) 



is the angular 



momentum (vector) operator for the i atom. Carte- 
sian components F^, Fy and Fz follow from this in an 
obvious way. Under physical conditions that preserve 
permutation symmetry of the label i, the collective spin 
of an initially polarized atomic sample can be restricted 
|l5j | to its maximum angular momentum shell. The 
associated Hilbert sub-space is spanned by eigenstates 
\F,M) of the collective spin observable F that satisfy 
F'^\F, M) = h'^F{F + l)\F, M), where F = Nf for atoms 
with individual spin /. 

It is natural to conceptualize the quantum state of 
such a system as a Bloch-like magnetization vector F = 
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FIG. 1: (A) Graphical representation of the spin-polarized 
atomic sample as a classical magnetization vector with trans- 
verse quantum uncertainty. (B) Schematic of the transverse 
quantum uncertainties for coherent and squeezed spin states. 
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FIG. 2: Schematic of an experimental apparatus for continu- 
ous measurement of collective spin in an Alkali atom sample 
based on polarimetric detection of a forward scattered probe 
laser. Information gained from the measurement can be used 
to achieve conditional spin squeezing. 

[{Fx), (Fy), {Fz)] plus a transverse uncertainty AFj^ as- 
sociated with the variances of Fx, Fy and F^ (see Fig.^. 
The transverse uncertainty AFj_ can never vanish since 
Fx, Fy and Fz do not commute; this constraint can be in- 
terpreted to mean that we can never have perfect knowl- 
edge of the orientation of the collective magnetization. 

Conditional spin-squeezing experiments operate on the 
principle that weakly measuring F^ gradually reduces 
its uncertainty below that of the initially-prepared spin 
state. A typical apparatus for such procedures is depicted 
in Fig. 121 Continuous measurement of F^ is implemented 
by passing a linearly polarized probe laser through an 
atomic sample prepared in an initial (approximate) co- 
herent state d6] by optical pumping |T3|. Qualitatively 
speaking, the atoms rotate (via state-dependent optical 
activity) the probe polarization by an amount propor- 
tional to Fz |l8|. In a quantum analysis the atoms and 
optical probe field evolve into an entangled state |3 ^ 
the result of this scattering interaction. Detection of the 
scattered probe field then yields information about F^ 
via these correlations P0II21II2I l2ll23| . 

The interpretation of existing squeezing experiments 
has (at least implicitly) assumed that polarimetric detec- 
tion of the forward scattered probe laser yields a detector 
photocurrent, 

ytdt = VMFzdt + dWt. (2) 

Here M is a constant (known as the measurement 
strength) that describes the rate at which photodetec- 
tion provides information about F^. The dWt are noise 
increments which exhibit Gaussian white noise statistics 
M[dWt] = and dW^ = dt 

The measurement strength, M, is the key parame- 
ter for predicting the degree of squeezing that can be 
achieved as a result of the measurement. It is thus 
important to determine Af in terms of physical prop- 
erties of the atomic sample and optical probe. While 
the form of Eg. ll^ has been derived in previous analy- 
ses P0II21I [2112a l2ll25l|. this has generally been done 
relative to a simplifying assumption that the atoms be- 
have qualitatively as spin-^ particles. But measurement- 
induced spin squeezing experiments have utilized Alkali 
atoms with higher spin [5|, |3, l26l | , and recent data show 



that the deviation from spin-i behavior can be signifi- 
cant |23|. While nonlinearities in the atom-probe scat- 
tering process are not always bad (proposals for capi- 
talizing on these effects for quantum state tomography 
are being explored 28] ), they do raise complications for 
spin-squeezing experiments by invalidating the form of 
Eq. ©. 

We find that the photocurrent in Eq. Q can be re- 
covered even for higher-spin atoms by suppressing ten- 
sor scattering interactions via a properly chosen exper- 
imental geometry. Using standard techniques ^29, S^l 
to address the atom-probe scattering physics, combined 
with a semiclassical treatment of the atomic magnetiza- 
tion vector, we derive an expression for the measurement 
strength, M, in terms of characteristic experimental pa- 
rameters. This allows us to obtain an expression for the 
photocurrent in terms of the duration of the measure- 
ment and the properties of the atomic system and the 
probe laser. We observe close agreement between our 
scattering model and data obtained using an apparatus 
of the type in Fig. 13 

Finally, we derive an expression for the signal-to-noise 
ratio of the measurement photocurrent which can be used 
to calculate the rate of spin-squeezing in experiments of 
the form shown in Fig. [3 The results we obtain are 
valid in the short-measurement limit in which atomic de- 
coherence due to scattering probe light in unobserved 
(non-paraxial) electromagnetic field modes [sj can be 
safely ignored. Current spin-squeezing experiments all 
fall into this regime where the degree of quantum un- 
certainty reduction is small compared to the Heisenberg 
spin-squeezing limit j^. 



II. CONTINUOUS MEASUREMENT AND THE 
MULTILEVEL ATOM-PROBE HAMILTONIAN 

We begin by considering the experiment depicted in 
the schematic of Fig.|5| An off-resonant linearly polarized 
probe beam is sent through a cloud of cold spin-polarized 
atoms. The forward scattered polarization state of the 
light is then detected using a polarimeter, consisting of 
polarization shifting waveplates, a polarizing beam split- 
ter, and two detectors.. 

In general, it is a rather complicated problem to pre- 
dict the output polarization state of the probe beam after 
it has passed through the spatially extended atom cloud. 
We can simplify the problem to one dimension by assum- 
ing that the beam is predominantly forward scattered 
due to the coherent re-radiation from a large number of 
atoms. This approximation can be extracted from a full 
three-dimensional model of the diffraction as considered 
in references |32, 33, 34] . Under this assumption, we 
only consider paraxial modes of the beam. Neglecting 
non-paraxial modes prevents us from computing the de- 
coherence rate of the atomic magnetization, but it does 
not limit our ability to analyze the dynamics in the small- 
decoherence (short measurement time) limit. 
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Even in the one-dimensional problem, the depth of the 
atomic cloud along the probe direction introduces further 
complications. To simplify further, we assume that the 
overall optical density of the cloud is small enough that 
the total rotation of the optical polarization state due to 
the atoms is small. This allows us to neglect propagation 
effects by which the atoms in the front edge of the cloud 
would see a substantially different input state than the 
back edge of the cloud. These effects lead to complicated 
multi-mode dynamics which are considered (along with 
the tensor polarizability) in 26]. 

Under these approximations, we approach the simpli- 
fied scattering problem as follows. The probe beam con- 
sists of two orthogonal polarizations and, for each polar- 
ization component, the continuous beam is divided into a 
series of distinct spatial traveling-wave modes, each with 
a length equal to the depth of the atomic cloud, L. Thus 
each atom interacts with a pair of orthogonal polariza- 
tion modes with the same spatial profile for a length of 
time St — L/c. This approach clearly avoids propagation 
effects by allowing all of the atoms to interact with the 
same modes simultaneously. After the two polarization 
modes corresponding to one time-slice have interacted 
with the cloud for the discrete time 6t , those two modes 
are detected with the polarimeter, a new time-slice begins 
to interact with the cloud, and the process is repeated, 
leading to a continuous measurement. More detailed ap- 
proaches to continuous measurement with discrete modes 
can be found in references p^ . l35l | . 

Now we define the electric-field and polarization oper- 
ators associated with each of these optical modes before 
considering the Hamiltonian interaction of probe beam 
with a single multilevel atom. This procedure is discussed 
with more detail in Appendix 1X1 



A. Probe Field Polarization States 

For each traveling-wave spatial mode, we consider the 
field operators 



and 



a_e_ + (2,6, 



(3) 



(4) 



where a[_ and a_ are Heisenberg-picture creation and 
annihilation operators for the z-axis propagating mode 
with left circular polarization and a^j^ and a_^_ are the 
creation and annihilation operators for right circular 
polarization. Each field operator implicitly refers to a 
single traveling mode as discussed above, and we neglect 
to provide indices for the modes unless they are required 
for clarity. The coefficient g — LLjo/{2eoV) is a form 
factor, V will be taken to be the volume of the atomic 
cloud, and e_ and e_|_ are the (complex) spherical basis 



vectors for left and right helicity. 

In the expansion of the polarizability Hamiltonian we 
get terms which can be recast as Schwinger boson oper- 
ators 

^0 ^ \ [ala^ + alaj^ (5) 
S'„ = - ( a\a_ + oLa, 



2 ('^y^y ^x^x) 



S,, — - { a^_a, — a\a_ 



I 
2 

= 2K'V-<'"x' 

These operators obey the usual angular momentum com- 
mutation relations and the components form a basis for 
the Stokes vector which is used to represent the polariza- 
tion state of the light. The quantity is proportional to 
the number of photons interacting with the atomic sys- 
tem in one time increment. On any given measurement, 
the quantity S'q and a single component of the Stokes 
vector representing the polarization state (e.g., S^) can 
be measured with an appropriate selection of polarization 
rotating waveplates situated after the atoms and prior to 
the polarizing beam-splitter. In the usual configuration 
(of Fig. (121)), S'x is measured without any waveplates, 
is measured with a half-waveplate that rotates the lin- 
ear polarization by 45-degrees, and is measured by 
adding a quarter-waveplate that completely circularizes 
linear polarized light. 

In the case where a full quantum mechanical descrip- 
tion is used, this choice of basis will change the nature of 
the information gained from the measurement which is 
then used to update the conditional collective quantum 
state describing the atoms. In other words, the choice of 
basis will lead to a different unravelling of the conditional 
dynamics. 



B. Scattering Hamiltonian 

We now introduce the polarizability Hamiltonian that 
determines the joint evolution of the single-atom spin 
and the polarization of the traveling-wave optical mode. 
Subsequently, we summarize the results from Appendix 
^ where we derive a more convenient and intuitive way 
of representing the irreducible components of the Hamil- 
tonian in terms of atomic spin operators instead of dipole 
operators. 

For a field which is off-resonant to the transition of 
interest, the usual dipole Hamiltonian can be approxi- 
mated and recast into a polarizability form. This can be 
derived, for example, by using adiabatic elimination un- 
der the assumption that the off-resonant field only weakly 
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populates the excited states. The polarizability Hamil- 
tonian 0,113, HI] is then expressed as 



H 



PfdPf'd^Pf 



hA 



(6) 



where we omit indices identifying the particular atom 
and spatial optical mode being considered. This defini- 
tion consists of several terms which are also defined in 
Appendix^ The negative and positive frequency probe 
field operators, E^^^ and E'^+\ describe the creation and 
annihilation of photons in the contributing probe modes. 
The atomic operators and d are the vector dipole 
raising and lowering operators. The ground and excited 
state angular momentum numbers are give by / and /' 
respectively. The probe detuning, Afj> = to — ^^fj', 
is defined as the difference between the probe frequency 
uj and a particular atomic resonance frequency. For the 
purposes of this paper we consider all of the population to 
remain in one ground state manifold (/ = 4 for Cesium, 
ignoring / = 3) and sum only over the excited states 
(/' = 2, 3, 4, 5). The operators P/ and Pf are projectors 
onto the ground and excited states respectively. 

This Hamiltonian has a satisfying physical interpreta- 
tion as a scattering interaction: the atom is first brought 
from its ground state to a virtual excited state via the 
raising operator, d^, by annihilating a photon from the 
probe field through Then, the temporarily excited 

atom returns to a (potentially different) ground state by 
emitting a photon into a (potentially different) scattered 
probe mode via d and E^^^. 

The central operator in the scattering Hamiltonian, 



(7) 



commonly called the atomic polarizability tensor, is a 
dyad involving vector operators |29j |. Thus a/j' is a 
rank-2 spherical tensor that can be decomposed into ir- 
reducible components. 
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The scattering Hamiltonian similarly decomposes into ir- 
reducible spherical tensor operators, 



where 



ijW=^E(- 

/,/' 



a 



fJ' 



hA 



e(+). 



/,/' 



(9) 



(10) 



The is a scalar contribution, H^^'> transforms as a 
vector, and //^^^ transforms as a rank-2 symmetric tensor 
in the group representation theory of SO {3). Were the 
atomic system composed of spin-i particles, it would be 
possible to neglect the rank-2 Hamiltonian [s^l (as will 
become explicit), however, we can not do so for higher- 
spin Alkah atoms [Ulllli^- 



The full Hamiltonian for the collective atomic spin re- 
sulting from N atoms is obtained by taking the symmet- 
ric sum of these single particle operators. 



C. Hamiltonian Decomposition 

Now we recast the single atom Hamiltonian of Eqs. 
(PI I10|l into irreducible terms involving only atomic spin 
operators fi and probe polarization operators St then dis- 
cuss each component in physical terms. The derivation 
of these expressions is sketched in Appendix IXI 



1. The Scalar Hamiltonian 

The scalar scattering Hamiltonian, i/^"^ , can be repre- 
sented as a product of operators on the separate atomic 
and probe field Hilbert spaces. This is accomplished 
by combining the expressions for the field mode oper- 
ators, Eqs. (O and with the rank-0 irreducible com- 
ponent of the atomic polarizability tensor. Evaluating 
this Hamiltonian using the form of the rank-0 atomic 
polarizability derived in Appendix 1X1 leads to the scalar 
scattering Hamiltonian, 



/' 



(11) 



where the constants ct^Pfi, defined in equation HA28() of 
the appendix, are related to the transition dipole matrix 
elements for the atomic hyperfine transitions. 

This rank-0 Hamiltonian couples the atomic identity 
operator 1^ to the field mode number operator and 
can be interpreted as an atomic state-independent light 
shift. It therefore affects both polarization modes of the 
probe field in an equivalent manner and will not infiuence 
the measurement process since it does not provide any 
state-dependent information. However, this Hamiltonian 
would be important if the measurement was meant to 
distinguish between populations across hyperfine states 
(e.g. / = 3 and / = 4 using homodyne detection) instead 
of across the sub-level populations within one hyperfine 
state (using polarimetry, as discussed here). This term is 
also of importance if the Hamiltonian is being considered 
as a spatially dependent potential for the atoms (e.g., in 
an optical lattice). 



2. Ttie Vector Hamiltonian 

The vector contribution to the atom-probe scattering 
Hamiltonian, can be evaluated in a similar manner using 
expressions for the rank-1 polarizability derived in the 
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appendix, 

H''' - 9T.ir^sj\- (12) 

Here, the vector polarizability constant, ctj^j-/ , is given by 

Eq. (|A29p , is the z-component of the (single-particle) 
atomic spin angular momentum. 

The rank-1 Hamiltonian can be interpreted as causing 
a differential phase shift on the two circular polariza- 
tion modes by an amount that is proportional to the 
z-coniponent of the atomic angular momentum. Thus 
the vector Hamiltonian leads to optical activity in the 
atomic sample and produces the familiar Faraday rota- 
tion effect often used to address continuous measurement 
of collective spin [H ^ [H IH IH . 




FIG. 3: Definition of the spherical coordinate angles used to 
describe the orientation of the collective atomic magnetization 
vector, F, relative to the fixed laboratory cartesian coordinate 
system. The polarization vector of the input probe light re- 
sides in the a;y-plane and forms an angle, (jip, with respect to 
the laboratory x-axis. 



3. The Tensor Hamiltonian 

Finally the tensor Hamiltonian, can be evaluated using 
expressions for the rank-2 polarizability derived in the 
appendix to give, 

+ Sy{fJy+fyL) 

+ 4(3/2 _/(/+i)i^)/3). 

(2) 

Here, the tensor polarizability constant, ctffi, is given 
by Eq. (|X30|) . 

The rank-2 Hamiltonian couples spin coordinates to 
the elliptical components of the probe laser field and 
produces a second-order light shift proportional to the 
atomic quadrupole moment. These terms vanish for 
/ = 1/2 (as can be seen by evaluating the operators 
within the parentheses above) but are non-zero for any 
higher spin number. For a linearly polarized input beam, 
the tensor term leads to an elliptically polarized scattered 
probe field [2^ |29|. The rank-2 interaction potentially 
limits the validity of any analysis of a continuous mea- 
surement of collective atomic spin in Alkali atoms based 
on the qualitative behavior of spin-i particles. 

D. Semiclassical Evolution of Probe State 

We can greatly simplify the dynamics by eliminating 
atomic evolution due to the probe beam and only consid- 
ering the evolution of the probe beam due to the atomic 
state. Under this semiclassical approximation, we replace 
all atomic operators with their expectation values with 
respect to an assumed fixed spin state. (This is the op- 
posite of the semiclassical situation often considered in 



atom-light interactions where the atomic system is con- 
sidered quantum mechanically while the optical beam is 
made classical.) For a large ensemble of atoms and small 
interaction times, fixing the atomic state will accurately 
reproduce the mean behavior of the measured photocur- 
rent corresponding to one of the Stokes vector compo- 
nents. This is confirmed experimentally in the next sec- 
tion, where the atomic state is fixed and adiabatically 
positioned with a magnetic holding field. The holding 
field serves to both position the atomic state and protect 
it from the infiuence of the probe light, such that the 
analysis of this section remains valid even for long inter- 
action times or large optical depth clouds. Ultimately, 
however, probe induced decoherence will dominate all in- 
teractions. In the final section, we then reconsider the full 
analysis including the atomic quantum noise (related to 
spin squeezing) for a particular alignment of the collec- 
tive spin state. 

We approximate the iV-atom Hamiltonian, Hni by re- 
placing the single-atom operators with their expectation 
value taken with respect to an optically pumped spin 
pointing with direction 9 and ip given in spherical coor- 
dinates (Fig. |3J). In other words, for an individual atom 
operator O/, we take 

6/^(O/) = (*(0,^)|d/|vi/(0,^)) (14) 

where |*(6',(^)) =exp[-if^(f]ex.p[-ify9]\f,f)z. 

The relevant operators from the Hamiltonian decom- 
position are given by 

(/J = fcosB (15) 
{fl -fy) = /(/ - 1/2) sin^ e cos 2^ (16) 
iUy-fyL) = /(/-l/2)sin2 0sin2^ (17) 

Within the semiclassical approximation, we obtain an ef- 
fective scattering interaction Hamiltonian that only in- 
volves operators on the probe field Hilbert space. Ignor- 
ing all terms proportional to Sq (because it commutes 
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H 



(1) 

N 



&(2) 



H 



N 



with each term of the semiclassical Hamiltonian) we have The quantity, (Jq , is the resonant atomic scattering cross 

section and A = 7rr^ is the cross-sectional area of the 
atomic sample. 

In Appendix IbI. the equations for a general rotation of 
(18) S about 7 are given. Here we specialize to the case where 
the input beam is linearly polarized in the x-direction 
such that (Sy) — (Sz) = 0. The output expectation 
values are then given by 



hxS^ + lySy+JzS^)-^ 



which leads to a rotation of the Stokes vector S about a 
vector 7 = [7x, 7j/, 7z] according to the evolution operator 



St 



exp 
exp 



(19) 



where dt = L/c is the interaction (transit) time of the dis- 
crete spatial modes of the probe beam across the atomic 
cloud of length L. The rotation vector 7 is defined by 



{S'J = (5,) (cos7 + ^(l-cos7) 
(SL) = (4)(-^sin7+^(l-cos7) 



(1 — cos 7) 



(26) 



{S'z) = (5.)psin7+^ 
V 7 7 



Taking the total rotation angle small (7^1) this be- 
comes (to second order in 7) 



- 7o/(/-l/2)sin2 0cos2vp^ 



(2) 



- 7o/(/-l/2)sin2 0sin2^^ 



= 70/ cos 6*^ 



-f' "oA/j- 



(20) 



(21) 



(22) 



Here we have normalized by the state-independent polar- 
izability constant (see Appendix 



such that aj y-z/ao is dimensionless. 
strength is represented by 



87r2 
(j||d||/> 



(23) 



(2j + 1) 
(2j' + 1) 



The rotation 



(5,) (1 - 7^/2 - 7,V2) 



('5'x) (7^ 



2 /2 _ „,2 

7y7 
2 

Izlx 



(27) 
(28) 

(29) 



In this semiclassical approximation, we have com- 
pletely neglected any evolution of the atomic state due 
to the probe beam. We demonstrate in the next section 
that the above model agrees well with experimental data 
when the spin state is fixed with a magnetic holding field. 



III. EXPERIMENTAL RESULTS 

In this section, we show that the model described above 
is consistent with representative data collected from our 
experiment with laser cooled Cs atoms and balanced po- 
larimetric detection of a forward-scattered, off-resonant 
probe laser field. 
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NgStao 



where we have used the field coeflacient g = LUo/{2eoV), 
the atomic resonance frequency loq, and the interaction 
volume (the volume of the atomic sample) V = AL. 

From an experimental standpoint, it is useful to note 
that 7o is directly related to the on-resonance optical 
depth OD of the atomic sample and the decay rate T 



via, 



where 



70 = ( 4 ) OD 



OD = TV 



A' 



CTO 



3A| 
27r ■ 



(24) 



(25) 



A. Experimental Apparatus 

Figure 01 provides a schematic of the major compo- 
nents of the experimental apparatus. Our single-particle 
Alkali atom spin system is the 62Si/2(f=4) ground state 
hyperfine manifold in ^^'^Cs with 4h of intrinsic angular 
momentum due to a combination of the i = 7/2 nuclear 
spin and the s = 1 /2 spin of an unpaired 6s valence elec- 
tron. We obtain cold atom samples from a 10~^ Torr 
background Cs vapor using standard laser cooling and 
trapping techniques by collecting more than 10^ atoms 
in a magneto-optic trap (MOT). Trapping beams are de- 
rived from a 150 mW injection-locked diode laser tuned 
(11-15) MHz red of the Cs 62Si/2(f=4)^62P3/2(f'=5) 
cycling transition. Each 35 mW trapping beam has 
an approximately constant intensity profile and a 2.5 
cm diameter. A 10 mW repump laser tuned to the 
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FIG. 4: Schematic of our experimental apparatus in which collective spin angular momentum of a cloud of laser cooled Cs 
atoms is measured by polarimetric detection of a scattered off-resonant probe laser. Ambient magnetic field fiuctuations are 
supressed by magnetic shielding and can be monitored with a fluxgate magnetometer (FG) situated nearby the atomic sample. 
Components not shown include the optical pumping laser (aligned along the laboratory a;-axis) and external trim coils used to 
zero ambient magnetic fields and their first order gradients. 



6^Si/2(f=3)^6'^P3/2(f'=4) transition is used to prevent 
atomic population from decaying out of the trapping cy- 
cling transition. 

Following the atom collection phase, the sample is sub- 
Doppler cooled '38] to a temperature of T ~ 10 iiK and 
the initial x-polarized spin state is prepared with a cir- 
cularly polarized 100 fiW optical pumping beam (pulsed 
for 2-4 ms) propagating along the a;-axis and tuned to 
the (f=4)— >(f'=4) hyperfine transition. A 100 mG mag- 
netic holding field is applied along the laboratory cc-axis 
to define the optical pumping direction. 

Continuous measurement of the polarized atomic en- 
semble is implemented with a nearly quantum shotnoise- 
limited probe laser that can be detuned from the 
6^Si/2(f=4)^6^P3/2(f'=5) Cs transition over a range 
A = ±1.4 GHz. The probe beam is linearly polarized by 
a high extinction Glan-Thompson prism prior to pass- 
ing through the cold atom cloud, and the orientation of 
the linear polarization vector with respect to the labora- 
tory coordinate system may be rotated via an input half- 
waveplate. The scattered probe field is detected with a 
polarimeter constructed from a Glan-Thompson polariz- 
ing beam splitter and a DC-balanced photodetector with 
>1 MHz measurement bandwidth. 

A computer controls the experiment timing and 
records the polarimeter output as well as diagnostic in- 
formation including background magnetic field fluctu- 
ations (measured with a flux-gate magnetometer) and 
atom number (measured by fluorescence imaging). The 
computer enables/disables the measurement by control- 
ling a shutter on the probe laser, constructed from a 
switched acousto-optic modulator, with 100 ns resolu- 
tion. Magnetic fields with magnitudes up to ~ 0.5 G can 
be applied in arbitrary (time-dependent) directions by 
driving 3 pairs of computer-controlled Helmholtz coils, 
oriented along the laboratory x-, y-, and z-axes, with a 
bandwidth of 1 MHz. 

Background magnetic field fiuctuations are suppressed 
through a combination of passive /i-metal shielding and 



field cancellation via external trim coils. Each atom 
preparation (trapping, cooling and optical pumping) and 
measurement cycle is synchronized with respect to the 
60-Hz building power lines to suppress the effects of in- 
duced magnetic fields. Slow magnetic drift due to natural 
and anthropogenic sources are cancelled by adjusting the 
external trim coils based on the output of the fluxgate 
magnetometer . 

B. Verification of the Probe Scattering Model 

Our model of the scattered probe polarization as a 
function of the orientation of the atomic magnetiza- 
tion vector was compared against experiment by observ- 
ing the polarimeter photocurrent as the orientation of 
the atomic polarization was varied according to differ- 
ent specifled paths in the laboratory coordinate system. 
This was accomplished as follows. An a:-polarized cold 
atom sample was prepared according to the description 
above and an a;-axis magnetic holding field of 100 mG 
was applied. At this point, the probe shutter was opened 
and the balanced polarimeter photocurrent was moni- 
tored while the orientation of the magnetic holding field 
was varied according to the specified path. The rate of 
change of the holding field orientation was chosen to be 
slow (ms) compared to the atomic Larmor precession fre- 
quency (hundreds of kHz) such that the atomic magne- 
tization vector adiabatically followed the path traced by 
the holding field. Furthermore the holding field was large 
enough to dominate the probe light induced dynamics at 
short times, but not so strong as to shift the levels sig- 
nificantly compared to the detuning. 

With a strong enough holding field, the spin state (and 
hence the semiclassical rotation vector 7) will be fixed 
across the spatially extended cloud. Because rotations of 
the Stokes vector about the same vector will commute, 
the semiclassical analysis of the previous section will be 
valid for even large optical depth samples where the total 
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FIG. 5: Comparison of our model of continuous measurement with photocurrents obtained from from the experiment with 
= 1 X 10^ Cs atoms in an r = 4 mm spherical trap and a P = 10 probe field blue-detuned from the (f=4)^(f' =5) D2 
hyperfine transition (Ao — 852 nm). Each trajectory is averaged 10 times. (A) For an input probe beam with x-polarization 
and a detuning of 150 MHz, Sy and Sz were measured for both the xz and xy trajectories (described in the text) resulting in the 
solid curves. All trajectory times are r = 2 ms, during which we observe some atomic decoherence which causes the prediction 
(dotted curves) to stray from the data. (B) As a function of probe detuning, we plot the peak of the Sy measurement (for the 
xz trajectory) which depends only on rank-1 terms and the peak of the Sz measurement (for the xy trajectory) which depends 
only on rank-2 terms. The predicted behavior (solid curves) shows good agreement with the data out to large detunings where 
the curves asymptote to the 1/A and 1/A^ lines provided to guide the eye. 



optical polarization rotation is significant. 

This process was performed for two different adiabatic 
paths on the atomic Bloch sphere: 

• a;z-Plane Rotation: the atomic magnetization 
follows a path beginning along the a:;-axis and ro- 
tates around the y-axis: 9 = n /2 ^ — 7r/2 with 
fixed = 0. 

• a;i/-Plane Rotation: the atomic magnetization 
follows a path beginning along the x-axis and ro- 
tates around the z-axis: — Q tt with fixed 
e = 7r/2. 

We chose these two trajectories because they highlight 
the different contributions from the rank-1 and rank-2 
scattering interactions, as seen from Eqs. (|20l 1211 I22|l . 
The ccz-plane trajectory, where = 0, virtually elim- 
inates the rank-2 tensor contribution to the photocur- 
rent leaving nearly ideal Faraday rotation. Conversely, 
the xy-plane rotation eliminates rank-1 contributions and 
produces elliptical scattered probe polarizations. 

1. Measuring the Scattering Probe Stokes Vector 

Fig.\^^ compares the measured polarimeter photocur- 
rents (solid curves) for these two adiabatic trajectories 
with those predicted by our atom-field scattering model 
(dotted curves). The input state was polarized in the x 
direction and because the total polarization rotation an- 
gle induced by the atoms 7 was small, we measured only 



the other two components Sy and with the appropri- 
ate arrangement of waveplates prior to the polarimeter. 
For measuring Sy a single half-waveplate is placed prior 
to the polarizing beamsplitter (PBS) to rotate the polar- 
ization by 45 degrees, and for S^ a quarter-waveplate is 
used to circularize the initial linearly polarized light. 

We now refer to Eq. and Eq. to explain the 
observed trajectories. For the xz trajectory, we have 
= such that Sy contains a large linear term in 7^ 

but S^ only contains terms quadratic in 7. Thus, for this 
path, the measurement of Sy leads to the top curve in 
Fig. which is proportional to the rank-1 polarizabil- 
ity, while the measurement of S'^, is much smaller and 
effectively zero. For the xy trajectory, we have jz — 
such that S^ contains a large linear term in 7j, but Sy only 
contains terms quadratic in 7. Thus, for this path, the 
measurement of S^ leads to the middle curve in Fig. 
which is proportional to the rank-2 polarizability, while 
the measurement of Sy is much smaller and effectively 
zero. The doubling of frequency between the two domi- 
nant curves is a direct consequence of the tensor nature 
of the rank-2 term. 

Note that there is some structure expected in the two 
curves (quadratic in 7) which are approximately zero, 
but these are more polluted by technical noise and do 
not reveal any essentially new information about the in- 
teraction. For the larger curves (linear in 7), deviations 
of the measured photocurrents relative to the predicted 
values seen in Fig. result mainly from the fact that 
the model does not consider the probe-induced damping. 
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The predicted curves use values for the atom number, 
trap volume, probe power and detuning consistent with 
independent characterizations of those parameters. The 
atom number and trap volume were obtained from fluo- 
rescence detection of the MOT and a CCD image of the 
atom cloud, and the resulting values, = 1 x 10^ and 
r = 4 mm, correspond to an optical depth, OD ~ 7, 
which is consistent with absorption measurements that 
we performed. Given our uncertainty in measuring the 
number of atoms, it can be inferred that our optical 
pumping efficiency in these (relatively) optically thin 
atomic samples is no less than 85% (but is more likely 
>90%) [3. 



2. Relative Scaling of the Scattering Terms with Probe 
Detuning 

As further verification of our scattering model, we in- 
vestigated the scaling of the rank-1 and rank-2 contri- 
butions to the polarimeter photocurrent as a function of 
the probe detuning. With reference to Fig. (El^), the 
magnitude of the vector and tensor scattering interac- 
tions were measured from the peak amplitude of the Sy 
measurement (for the xz plane rotation) and the ampli- 
tude of the 5*2 measurement (for the ccy-plane rotation) 
respectively. This plot compares these measured signal 
amplitudes (stars) with those predicted by our scatter- 
ing model (solid curves) for detunings (with respect to 
the (f=4)^(f'=5) hyperfine transition) ranging from 150 
MHz to 1.05 GHz. 

The fact that multiple excited state hyperfine levels 
participate in the scattering interaction is evident from 
scalings which are not constant in A^^ or A^^. As sup- 
ported by our full model of the scattering interaction, we 
observe no qualitative difference in the continuous mea- 
surement for probe detunings smaller than the hyperfine 
splittings. This suggests that conditional spin-squeezing 
experiments can be performed with small detunings pro- 
vided that the probe intensity is weak enough that the 
small decoherence requirement is satisfied. 



IV. SPIN-SQUEEZING WITH MULTILEVEL 
ATOMS 



Until this point we have considered only the semiclas- 
sical evolution of the optical probe beam due to an en- 
semble of atoms with a fixed atomic spin state. Now 
we consider a different experimental scenario appropri- 
ate for preparing conditional spin-squeezed states of the 
atomic ensemble. As opposed to the previous situation, 
we remove the adiabatic holding field which makes spin- 
squeezing impossible as it will cause undesired mixing of 
the squeezed and anti-squeezed components perpendic- 
ular to the mean spin. Although the holding field may 
serve to validate the previous semiclassical analysis for 
longer times by eliminating the probe-induced evolution 



of the atomic state, this analysis is still valid for small 
times and weak interactions without a holding field. Thus 
our goal is to derive the small time signal-to-noise ratio 
by deriving the signal strength from the previous section 
and comparing this to the optical shotnoise. We then use 
the signal-to-noise ratio to predict the rate of squeezing 
in a typical experimental configuration where the tensor 
terms can be ignored. 

Considering only the relatively low optical density 
limit, the measurement of Sy will result in Eq. H28|) . Now 
we wish to re-write this equation in the form of Eq. (j^J 
including measurement noise. It is readily shown that 
all terms not linear in Fz vanish in Eq. (|28|l provided 
that 9 = tt/2 and = 0. That is, a pure Faraday ro- 
tation Hamiltonian is recovered when the atomic mag- 
netization vector is oriented along the x-axis. However, 
rotating F in the xy-plane results in elliptically polar- 
ized scattered probe light, and moving out of this plane 
results in nonlinear atomic dephasing due to scattering 
terms which are quadratic in the single-particle spin op- 
erators, fz- These adverse effects are avoided for the 
experimental geometry where F is collinear with the x- 
axis. Fortunately, spin-squeezing experiments are easily 
operated under such conditions j7(. 

Taking the input probe field to be in an x-polarized 
optical coherent state, and considering the small 7 limit. 
Equation (|28|l leads to a semiclassical photocurrent (with 
units of optical power) of the form. 



yt 



vVSFz + ^Ct 



(30) 



where we have made the substitution, hNf cos 9 
(refer to Eq. (|15|l 'l. and included the photodetector quan- 
tum efficiency, 77. Note that we have introduced Ct which 
represents optical shotnoise. We have also introduced a 
constant, S, the scattering strength, 



IpCTO [ J 



(31) 



that depends up the probe intensity, Ip — P/A, deter- 
mined by the coherent state amplitude, P = 2huj\f3\'^ and 
cross-sectional area, A — 7rr^ (for a mode-matched probe 
laser). It is useful to note that the scattering strength 
has units of W^/fi^ (power squared per ti^) and charac- 
terizes the degree of coupling between the atoms and the 
probe field; ^/S quantifies the polarimeter optical power 
imbalance per unit spin (as F^ has units of H). 

Our expressions are similar to previous results 1^ 123 . 
j3.9| in that it appears as a Faraday rotation signal. How- 
ever, our specific expressions for 7x, 7y and 7z account 
for the detailed hyperfine structure of the atomic excited 
states, including the fact that the oscillator strengths and 
signs of the contributions from different participating ex- 
cited states are not equal, and doing so is required for 
quantitative agreement between theory and experiment. 

To arrive at an expression for the measurement 
strength, M, as defined in Eq. 0, we must consider the 
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variance, AC^, of the white noise increments For an 
optical coherent state |4^.41j, this noise variance is given 
by the famihar optical shotnoise expression, 



(32) 



which has units of W^/Hz (power squared per frequency). 
Comparing the semiclassical photocurrent of Eq. (|30|l to 
the photocurrent of Eq. (01, the measurement strength is 
seen to be given by the ratio 



M 



S 



1 



(33) 



AC 2h? 

where we have defined the reciprocal scattering time as 



4 ^ aoA/j' 



(34) 



which is essentially the rate that probe photons are scat- 
tered by the atomic system. This expression is similar to 
that derived in Ref. [l8l| . 

Now consider a measurement of by Eq. In the 
small time limit where probe induced decoherence can 
be neglected, the full quantum filter describing this mea- 
surement is equivalent a classical model in which Fj. is 
simply a random constant on every trial drawn from a 
distribution with variance equal to the quantum vari- 
ance of {AF^)q l4a|. Then the generally complicated 
full quantum filter |25| is equivalent to linear regression, 
or fitting a constant to the noisy measurement record in 
real time. In essence, the optimal filter serves to average 
away the optical shotnoise to reveal the underlying value 
of i^z. Under these statistical assumptions, at small times 
the quantum uncertainty is given by 



l + 77(AF2)oAfT- 



(35) 



This can be shown either with the full quantum filter 
or by using the equivalent classical model combined with 
Bayesian estimation (from which a Kalman filter or linear 
regression can be derived). 

These concepts are illustrated by the simulated mea- 
surement trajectory in Fig. The plot begins with the 
probe laser turned off, during which all necessary state 
preparation of the atomic system such as atom trapping, 
cooling and optical pumping into an x-polarized coher- 
ent spin state is performed. Once the probe light is en- 
abled at t = 0, the photocurrent acquires a mean offset, 
rj^/SFz, proportional to the spin measurement outcome, 
Fz, but this mean value is masked by photocurrent noise. 
At short times, the signal is overwhelmed by local statis- 
tical fluctuations; however, averaging the photocurrent 
suppresses the uncertainty in the mean signal by inte- 
grating away the white noise, illustrated by the dotted 
lines in Fig. |51 

If we define the signal to noise ratio as 




0.5 
Measurement Time (X, arb. units) 

FIG. 6: Simulated photocurrent (77 = 1) for a continuous 
measurement of atomic spin angular momentum via balanced 
polarimetry. At the onset of the measurement, t — 0, the 
photocurrent assumes a mean offset proportional to the z- 
component of the spin, but this offset is masked by white 
noise due, in part, to optical shotnoise on the probe laser. 
Filtering the photocurrent gradually reduces the uncertainty 
in the photocurrent offset and produces spin-squeezing. 



we can then express the degree of squeezing (ignoring 
decay of the F^) as 



W 



1 

1 + snr2 



(37) 



Using {AF^)o = fi^Nf /2, we can express the signal to 
noise ratio as 



SNR2 = tjOD- — 

4 r. 



(38) 



SNR2 = ri{AF^)oMT 



(36) 



To keep this expression valid we must have r <C Ts, so 
our only recourse to creating large amounts of squeezing 
in free space is to increase both the quantum efficiency rj 
and the optical depth OD as much as possible. 



V. CONCLUSION 

In this work, we have derived the most useful form 
of the polarizability Hamiltonian describing the realistic 
measurement of an ensemble of multilevel alkali atoms 
with an off-resonant probe beam. We then showed that 
this model was consistent with experimental observations 
in the semiclassical limit where the atomic state was adia- 
batically directed with a strong magnetic field. We found 
that an adequate comparison was only possible after in- 
cluding all relevant hyperfine transitions including their 
relative (non-unit) oscillator strengths in our model of 
the atomic physics. 

We then developed a model for describing conditional 
spin-squeezing in Alkali atoms. Detailed investigation 
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of the atom probe scattering physics indicates that it is 
possible to ehminate unwanted tensor components of the 
atomic polarizabihty by adopting a suitable atomic and 
optical polarization geometry. This includes the elimi- 
nation of dephasing due to the quadratic light shift p7| 
without sacrificing a fixed laboratory coordinate system 
for the measurement. Moreover, we found that condi- 
tional spin-squeezing experiments could be performed at 
small optical detunings without a substantial change in 
the form of the photocurrent or filtering approach. 
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APPENDIX A: THE IRREDUCIBLE 
REPRESENTATION OF THE POLARIZABILITY 
HAMILTONIAN 

In this appendix, we derive the irreducible components 
of the polarizability Hamiltonian, Eqs. (|11II12|[T^ . from 
the less useful form of Eq. ©. We begin by clarify- 
ing notation used for the spherical basis and the spin 
states of the Alkali atoms. Then we discuss properties of 
the polarizability Hamiltonian and the dipole operator 
before detailing the decomposition and re-formatting of 
the Hamiltonian into its irreducible form. 



1. Spherical basis 

The spherical basis is the preferred basis when dealing 
with atomic transitions due to its symmetry properties. 
The basis is defined by the transformation from Cartesian 
coordinates 

e+ = -{S., + iey)/V2 (Al) 
eo = 

Thus elements in the spherical basis have the properties 
e;=e_,(~l)« (A2) 

— * — r 

and for an arbitrary vector A we have = • A so 
that A = A,el = A^-g. 

2. Alkali spin states 

We represent the internal state of the atom in terms of 
the (Zeeman degenerate) atomic hyperfine states, |/, m). 
Here / and /' are the total spin quantum numbers for the 
ground and excited hyperfine levels respectively while m 
and m' are their projections on the z-axis. That is to 
say, I/, m) are eigenstates of the total atomic angular 
momentum, 

f = s (g) i/g, -I- (g) i «) li -h i^^ii (A3) 

where s, 1, and i are respectively the electron spin, orbital 
angular momentum, and the nuclear spin. The quantum 
numbers, /, and to, are defined in the usual manner, 

fV,m) = h^f(f + \)\f,m) (A4) 
fz\.f,m) = hm\f,m) 



We use the notation that are in the spherical basis 
/±=T(/.±*/y)/V2. (A5) 

It will also be useful to define a projector onto the ground 

state / 

P/=^|/,TO)(/,m| (A6) 

and a projector onto the excited state /' 

Pj,=Y,\r,m'){f,m'\. (A7) 

3. Hamiltonian approximation 

We begin with the single-particle dipole Hamiltonian 
H = — d • E. The dipole operator d — efe can be split 
into its raising and lowering components 

d = d(-)+d(+) (A8) 
/,/' 
/,/' 

and the electric field operator can be split into rotating 
and counter-rotating terms 



E = 


e(-) 


f e(+) 




e(-) = 








E(+) = 




a_e_ 4 





After using the rotating wave approximation and one of 
many available perturbation expansion techniques (e.g., 
adiabatic elimination) we arrive at the familiar polariz- 
ability Hamiltonian 0, [s^, ll^ , 

H = yE(-).-^.E(+) (AlO) 
SJ' 

where the atomic polarizability between a particular 
ground state (/) and excited state (/') is defined as 

"/,/' = ^'/dP/'dtP/ (All) 
= EEEl/'"^")(/'"^"|d|/'-"^') (A12) 

m m' m" 

x(/',m'|dt|/,m)(/,TO|. 

This expanded expression involves dipole operator ma- 
trix elements of the form, (/', m'\dq\f^ to) where |/, m) is 
a Zeeman sub-level in the ground-state hyperfine man- 
ifold, |/',m') is a virtual state in the excited hyperfine 
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manifold, and q = 0, ±1 labels the helicity of the electro- 
magnetic field. 

The above notation is complete, but for the rest of this 
appendix we work with only one particular /, /' combi- 
nation and remove the subscripts with the simplifying 
notation change 



afj, a 



(A13) 



However, when the complete Hamiltonian is considered, 
the summation over all possible /, /' combinations is re- 
established. 



4. Matrix element decomposition 

In order to work with the above expressions, it is ad- 
vantageous to simplify the dipole matrix elements as 
far as possible. By employing the Wigner-Eckart the- 
orem, the angular dependence of the matrix element, 
(/', TO'|d|/, m) can be factored into the product of a 
Clebsch-Gordan coefficient and a reduced matrix ele- 
ment, 

(/,m|d;|/',m') - (/,m|l,g;/',m-q>(/||d||/'). (A14) 

Since the dipole operator acts only on electronic degrees 
of freedom, it is further possible to factor out the nuclear 
spin degrees of freedom via the explicit coupling, 

(/||d||/'> = (_i)/'+^+»+V(2/' + l)(2j + l) 

^'|01|d,||/) (A15) 

where i is the nuclear spin quantum number, j and f 
are the ground and excited state fine structure quantum 
numbers, and do is the dipole operator with respect to 
the electronic degrees of freedom. 



5. Tensor decomposition 



From Eq. (17-89) of reference [43 we see that we can 
form an irreducible tensor, Zm j from a linear combina- 
tion of tensor operators tlq'^^ and vj," ■* via the definition 

= E q; q'\j, m) (A16) 

where {K,q; k' , q'\j,m) are Clebsch-Gordan coefficients. 
This expression can then be inverted using 



We now specialize to the case where z'}^ = Tm\ U ~ 
d, and V = d^. Because we are creating a dyad (with 
two vectors), we have n = k! — 1. Inserting these above 
gives the definition 



f^^^ =Y,d,d\,{l,q-l,q'\j, 
and the inverse 



(A18) 



(A19) 



We can use this latter expression to write the polariz- 
ability as 



a = 



(A20) 
(A21) 



q,q' 



Y,Y.^;e^q'^^\l,q;l,q'\j,m) (A22) 

J>™ 9,9' 

(A23) 



where 



(A24) 



Filling in these Clebsch-Gordan coefficients explicitly, 
we get 



a 



(0) _ rp{0) 



a 



(1) 



^(1) 



1 1 . 

V2 + V2 



T3 



Gne_ — —G_Gn 



(2) _ rp(2) 

— J. n 



71 " 71 

7!'^°'^" + 71 



-T 



(2) 



V2^° + ' V2 
x/2 " V2 



+fi'^ [e;e;] 

+^1^ [e*_e*_] 



(A25) 



Furthermore, using the definition of Tm^ and filling in 



the Clebsch-Gordan coefficients explicitly, we get 



T, 



(0) 



(1) 



d+di - d_di 



(A26) 
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fii) 



7^(2) _ 
J- n — 



fi(2) 
fi2) 



T 



T 



(2) 
+2 

r(2) 



(d:i'_ + 2dodJ + d_d 
\/6 ^ 



(do^"+ + d+do 



Note that several standard references (including refer- 
ences 43j 44J) contain an error in the prefactor of the 
j = term and in the sign of the j = 1 terms. However, 
the fundamental definitions of Tj,P and its inverse above 
are valid. 

Using recursion relations for the Clebsch-Gordan co- 
efficients we can recast the tensor operators in terms of 
more intuitive / operators [2^ |45| 



r(o) 



f(i) _ 

^+1 - 

fii) ^ 

^0 — 

fi2) ^ 

rp(2) _ 

^+2 - 

fi2) ^ 



-4"},i,/V3 

2 



(A27) 



~af},{3f^ -/(/ + l)i^.)/%/6 



'V2/+ (/. 



(2) 



i,/2 



(2) ?2 
(2) ?2 



Here we have defined 



(0) 



af ((2/ - 1)<5;^1, + (2/ + l)<5f 

+ (2/ + 3)<5^+i) (A28) 



a 



(1) 
/,/' 



2/ + 1 



>,(2) _ 



1 



(A29) 



(A30) 



These definitions have been chosen to make the each of 



the quantities 



E 



> 



/,/' 



(A31) 



for Afji ^ for each term j. 
polarizability constants 



We have defined the 



a 



f 



and 



i2j + 1)2 



87r2 
(j||d||/) 



1 J f 
» /' / 



(A32) 



(A33) 



^ (2j + 1) 
(2/ + 1) 



which involves the atomic the spontaneous emission rate, 
r, and transition wavelength, Xq. 

Now, to complete the derivation, insert Eqs. l|J27)l into 
the polarizability components of Eqs. (|A25I) . then insert 
this and the definition of the electric field, Eq. HA9|I . 
into the Hamiltonian, Eqs. H9I10|I . Expanding, using the 
properties of the spherical dot product, and the Stokes 
component definitions (Eqs. |SJ, and summing over the 
/', we get the final expressions used in the text (Eqs. 1111 



APPENDIX B: ARBITRARY VECTOR 
OPERATOR ROTATIONS 



Here we are interested in evaluating the general opera- 
tion of rotating a vector about an arbitrary direction by 
an arbitrary amount in order to determine the semiclas- 
sical evolution of the probe light as used in Eqs. (|25|l . 

Consider the rotation of the vector spin operator 



(Bl) 



in Cartesian coordinates about an arbitrary direction 
n = [7rr,7y,7z]/7 by the angle 7 = ^ll+l^+ll- This 
rotation can be represented in the Heisenberg picture as 

S[ = ilS^U^ (B2) 

where 

U = exp[-i7S' • n] = exp[-i(7j;S'j; -|- 7^6"^ + 7^6"^)] (B3) 

The S[ can be derived explicitly using the following equa- 
tion for the arbitrary rotation of any vector 

5- = {S ■i)cos"f + {n-i)[n- S){1 ~ cos"i) 

-I- ((n X ?) • ^) sin7 (B4) 
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Expanding and rearranging terms we get 
S', = (^^(l-cos7) + cos7 



Ixly 

y » ^2 



(1 — COS 7) H sin 7 



+5. (^(l-cos7)-^sin7 



\ 7 



(1 — cos 7) — — sin 7 



+Sy ( -| (1 - cos 7) + cos 7 



(B5) 



+5, ('2^(l-cos7) + ^sin7 



V 7 



(B6) 



5; = 2f^(l_cos7) + ^sin7 
+5j'^(l-cos7)~2:^sin7 



7 



+5"^ ( ^(l-cos7) + cos7 



(B7) 



These equations can be specialized to Eqs. (|2mi which 
describes the experimental situation considered in this 
work. 



